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Abstract 

The paper concerns upper and lower estimates for the number of negative eigen- 
values of one- and two-dimensional Schrodinger operators and more general oper- 
ators with the spectral dimensions d < 2. The classical Cwikel-Lieb-Rosenblum 
(CLR) upper estimates require the corresponding Markov process to be transient, 
and therefore the dimension to be greater than two. We obtain CLR estimates in 
low dimensions by transforming the underlying recurrent process into a transient 
one using partial annihilation. As a result, the estimates for the number of negative 
eigenvalues are not translation invariant and contain Bargmann type terms. We 
show that a classical form of CLR estimates can not be valid for operators with re- 
current underlying Markov processes. We provide estimates from below which prove 
that the obtained results are sharp. Lieb-Thirring estimates for the low-dimensional 
Schrodinger operators are also studied. 

Key words: Schrodinger operator, negative eigenvalues, CLR estimates, Lieb-Thirring 
estimates. 
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1 Introduction 



Let Nq{V) = #{Aj < 0} be the number of non-positive eigenvalues of a Schrodinger 
operator 

H — —A - V(x), V>0, (1) 

on R d or Z d . Everywhere below we assume that the potential is non-negative. The 
standard approach to Cwikel-Lieb-Rosenblum (CLR) estimates for N (V) (see [6], [T5] - 
[18], [23], [22]) requires the Markov process x(t) which corresponds to the unperturbed 
operator H Q = —A to be transient (non-recurrent). The transience in the lattice case 
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means that the expectation of the total time the process x(t) spends in the initial point 
is finite. The latter is equivalent to the condition 



Po(t, x, x)dt < oo, (2) 
where po(t,x,y) is the fundamental solution of the corresponding parabolic problem 

-J? = &p , t > 0, p (0,x,y) = 8y(x). 

In the continuous case, one needs to talk about the time spent in a neighborhood of the 
initial point (not at the point itself), and the transience means that 



oo 



Po(t,x,y)dydt < oo, 

'o Jn 

where Q is a neighborhood of the point x. This condition implies that 

Po(t, x, x)dt < oo, (3) 

Conditions ([2]),([3]), obviously, do not depend on x. The integrals (T5]), ([3]) diverge for 
recurrent processes. The CLR estimates are valid for more general operators than (TjQ) (see 
[2U [25l [19] and references there) , but usually the transience is an essential requirement 
when these more general operators are considered. 

Recall one of the forms (not the most general) of the CLR estimate for Schrodinger- 
type operators H = H — V(x) on L 2 (X,B,fi) where X is a complete cr-compact metric 
space with Borel cr-algebra B(X) and a cr-finite measure fi(dx). Let H be a self-adjoint 
non-negative operator such that the operator — H is the generator of a Markov semigroup 
P t acting on C(X). Let po(t,x,y) be the kernel of P t , i.e., po is the transition density of 
the underlying Markov process x(t). If the process x(t) is transient, then 

N (V)<^— [ V(x) [ p (t,x,x)dtfi(dx), V>0, (4) 
c{a) J x J« 

V(x) 

where a > is arbitrary and c(a) = e~°~ J °° ze z ^ z ■ Another widespread form of this 
estimate 



N (V) <C d / V*(x)dx, V>0, (5) 

jR d 

for the Schrodinger operator in R d , d > 3, follows immediately from (Tj0) since po(t, x, x) = 
Cdt~ d l 2 in this case. This fact stresses the importance of the transiency requirement. It 
also explains the reason for imposing the assumption d > 3 in papers on CLR estimates. 

This paper concerns the CLR-type estimates for operator (DO) in dimensions d = 1 and 
2 where the Markov process is recurrent and estimates (H]),© fail. We also will consider 
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the fractional-dimensional case d < 2 presented by the Bessel operators and the fractional 
degrees of the one- dimensional discrete Laplacian. As usual, the Laplacian on the lattice 
Z d is defined as follows 

AV(a;)= (V^')-^))- (6) 

x':\x— x'\ = l 

The literature on the negative spectrum of Schrodinger operators is so extensive that 
we are going to mention here only some of the most closely related papers. The quazi- 
classical asymptotics for N (aV), 

N (aV) ~ c(d, V)a d/2 , d > 1, a -»■ oo, 

is valid [22J in dimensions cf > 3 when V G L d ^ 2 (R d ), and the CLR estimate (EJ is 
in agreement with the quazi-classical asymptotics. One could expect that N (aV) ~ 
c(V)a, a — > oo, when d = 2 and V G L 1 (i? 2 ). However, while this is true [22] for "good" 
potentials V, this asymptotics is not valid for arbitrary V G L l (R 2 ) [3] , [2] , [25] . One of 
our goals was to establish an estimate for Nq(olV) of the first or "almost" first order in 
a when d = 2. Thus we are not going to consider the Birman-Schwinger estimates pQ, 
[27], [22] which are not sharp when a ^> 1. The latter estimates, known for d > 3, were 
extended recently for low dimensions [31] . but they have order a 2 when d = 2. We will 
mostly focus on the two-dimensional case since one- dimensional problems are very specific 
and can be studied by a variety of methods. However, our general approach allows us also 
to obtain new one- dimensional results (for example, an improved Bargmann estimate and 
the Lieb-Thirring-type estimate with 7 < 1/2). 

Our two-dimensional estimates for jVo(aV) have order a, a — > 00, in the lattice case 
and order a In a for operators in L 2 (R 2 ). Simpler estimates of order a In a for N (aV) 
were obtained in [26], [21], [14] for two-dimensional operators with a central potential 
V(x) = V(\x\). Using a fundamental one-dimensional estimate of the Leib-Thirring sum 
from [TT] a remarkably simple formula is obtained in [14] : 

N (V)<1 + ^- [ V(\x\)\\n\x\\dx + ^-= [ V(\x\)dx. 
^ Jb 2 7tV3 Jr 2 

Unfortunately, this result is valid only for central potentials (a counter-example can be 
easily constructed if V is not central). The authors also formulated an elegant conjec- 
ture for general V. The conjectured estimate contains a term with a central decreasing 
rearrangement of the potential. An estimate of order a in two dimensions for general 
potentials was obtained in [30J in terms of the local Orlicz norms of the potential (an 
earlier exposition of the author's technique and its development can be found in [1], [5]). 
Our estimate in the continuous case is logarithmically weaker when a — > 00, but has a 
much simpler form. 

Our approach to treat the recurrent operators is particularly simple when the dimen- 
sion d is less than 2 and for the lattice two-dimensional operator. In these cases we replace 
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the unperturbed operator Hq = —A by its rank one perturbation Hi for which the killing 
(annihilation) of the corresponding Markov process in a single point is introduced. The 
idea to use the rank one perturbation in this problem is going to B. Simon [29J and was 
used in [31] . The difference is that we combine it with the CLR rather than with the 
Birman-Schwinger estimates, and this allows us to get the result in several lines. After 
the killing is imposed, the Markov process generated by Hi is transient. The standard 
CLR arguments can be applied to Hi — V, and therefore 

N (V) < / V(x) [ Pl (t,x,x)dtfi(dx) + 1, (7) 

where pi is the solution of the problem 

^ = -#iPi, t > 0, pi(0,x,y) = 5 y {x). 

Let us note that po and p\ are integrable near t — if d < 2 and in the lattice case. 
Thus one can use the estimate ((7j) with a = (and c(0) = 1). This implies the following 
particular version of ((7|): 



N (V) < / V(x) / p 1 (t,x,x)dtfj l (dx) + 1. (8) 
Jx Jo 

Consider, for example, a lattice Schrodinger operator in dimensions one or two. Then Hi 
is obtained by imposing the Dirichlet boundary condition at one point, for example, at 
x — 0. We will show that 



^2pi(t,x,x) = \x\, d=l, (9) 
^2 P i(t,x,x) < Cln(2 + |x|), d = 2. (10) 



Thus, OH]) implies 



A^o(^) < ^2\x\V(x) + l, d=l, (11) 



z 



N (V) < Cj2 ln ( 2 + \^\)V{x) + l, d = 2. (12) 

z 2 

The continuous analog of the estimate f lTT]) (for the operator in L 2 (R)) has the form 

iVo(^) < / \x\V{x)dx + l, (13) 
Jr 

and coincides with the well-known Bargmann estimate (see [22]). The Bargmann estimate 
has a wrong scaling order (a instead of y/ot). The estimates of order y/a can be found in 
see also the Calogero estimate [22] for monotone potentials. 
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Inequality (j7|) with a > leads to the following refined Bargmann's estimate for 
operator in L 2 (R) 

N {V) < -^—[ / \x\V(x)dx + — = / x 2 V 3/2 (x)dx] + 1, d = l, (14) 

c \?) Jx*V(x)>o V an Jx 2 V(x)<u 

with the some c(a) as in ([7]). Note that the Bargmann's estimate f|T3|) as well as other 
mentioned above one-dimensional estimates do not provide any information in the case 
of the potential 

V(x) = 0( - r ), |ac|-)-oo, d = l, 
x l In |x| 

(the integral in (TT3T) diverges), while the refined formula (fT4"|) shows that iV (V) < oo for 
this type of potentials. 

The most technically difficult part of the paper concerns the Schrodinger operator in 

R 2 : 

H = -A - V{x), H = -A, V>0, xe R 2 . 
A rank one perturbation approach does not work here. We replace H by 

Hi = —A + q(x), q = 1 for \x\ < 1, q = for \x\ > 1. 

The operator H = —A + q(x) — (V(x) + q(x)) can be considered as the perturbation of 
Hi by the potential V + q. The Birman-Schwinger principle implies that 

N {V) = N (V + q; Hi) < N (2V; Hi) + N (2q; H x ). 

where the second argument of the function A^ is the unperturbed operator. We chose q 
to be so small (q < 1) that No(2q; Hi) < 1, and therefore 

N (V)<N (2V;Hi) + l. 

We show that operator Hi is transient, and moreover the following non-trivial estimate 
holds: 

C(a)ln 2 (2 + Id) , , . . l2l . ., 

Pi{t, x, x) < 5 when t > ^{x) — max(l, a\x\ ln|x|), a > 0. 

tin t 



This implies that 



N (V) < -3— / V(x) I p 1 (t,x,x)dtdx + 1 
c(o-) JR2 



V(x) 



<Ci(a) [ Z^Ll n 2 (2+|x|)dx + C 2 (a) / Vln^^dx+1. (15) 
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A similar refined Bargmann-type estimate will be proved in the two-dimensional lattice 
case. 

Let us stress that all the above estimates are not translation-invariant unlike the case 
of Schrodinger operators in dimensions d > 3. The following arguments show that the 
estimates of the form 



for operators on Z d and similar estimates in the continuous case can not be valid in 
dimensions d = 1,2. Consider the operator H = —A — sS(x), x G Z d , d = 1 or 2. This 
operator with an arbitrary small e > has exactly one negative eigenvalue [H]. One can 
take an arbitrary sequence e = e n , n — 1, 2..., for which 'Y^ in s] l < 1 and choose a sequence 
x = x n to be so sparse that the eigenfunctions of the operators H n = —A — e n 8 Xn (x) 
are practically orthogonal, and the operator H = —A — 'Yl in £ n ^x n { x ) has infinitely many 
negative eigenvalues [33] . Then the left-hand side of ( |T6|) is infinity and the right-hand side 
does not exceed 1 + c. Similar arguments will be provided for Schrodinger operators on 
general lattices with the recurrent Markov process generated by the unperturbed operator 
(see section 6). 

We will also obtain the estimates on N (V) from below. The following results show 
that estimates (ITT]) - (ITS"]) are sharp (relatively to the decay of the potential at infinity). 
Consider the lattice Schrodinger operator on Z d . Then iV"o(V) = oo for two-dimensional 
operators with any potential V such that 



Similar results are valid in the continuous case. Note that a much stronger result is known 
[TO] for the continuous two-dimensional Schrodinger operator: 



Let us turn to non-integer dimensions. We will study two types of operators with 
fractional spectral dimension: Bessel operators and fractional powers of the lattice Lapla- 
cian. These are the operators on the half line R + and R, respectively, but their negative 
spectra behave as for multi-dimensional operators of dimension d that is not necessarily 
an integer. The Bessel operators are defined by 




(16) 




and for any one- dimensional operator when 




for some e > 0. 





d-1 d 



r dr 



■d-1 



1 d-1 



r 



(r 



d 
dr 



), d > 1 can be not-integer, 



(17) 
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in L 2 ([0, oo), r d ~ 1 dr). Operator Bd is selfadjoint when d > 2. If d is strictly greater than 
2, then the diffusion process bd(t),t > 0, generated by S^, is transient and the classical 
CLR estimates imply 

POO 

N (V) < C(d) / V d l 2 {r)r d - X dr. 
Jo 

If d = 2, then the process bd(t) is recurrent, and the estimates for Nq(V) can be derived 
from the results obtained in this paper for two-dimensional Schrodinger operators. 

The process b d {t) as well as the operator B d are not determined by (IT7|) if d < 2. 
One needs to add a boundary condition at r = 0. The operator B d and the process bd(t) 
become well-defined if the Dirichlet boundary condition is imposed at r = 0, and the 
following analog of (TH|) is valid in this case: 

N (V)< Cl (a) I V(r)r 2 ~ d dr + c 2 (a) f V 2 ~ d/2 (r)r 4 - 2d dr. 

Jr:r 2 V>(T Jr:r 2 V<a 

Somewhat similar results are valid for the fractional powers of the lattice operator. We 
will only consider the powers (— A) a of the one-dimensional lattice Laplacian. It will be 
shown that the Bargmann-type estimates are valid in this case with the constant d = 1 /a 
playing the role of the spectral dimension. Note that the Markov process generated by 
the operator — (— A) Q is not local (it has a positive probability of jumping to any point of 
the lattice). The results for the operator(— A) a will be obtained as a consequence of the 
results for more general Schrodinger operators on discrete graphs which will be studied 
in section 6. 

The final section of the paper is devoted to estimates on the Lieb-Thirring sums 

Sy(V) = £ \Xj\. 
i:\i<0 

The classical result [18J for the Schrodinger operators in R d has the form: 

S 7 (V)<c djJ f V^dx, ~ + 7>l- (18) 

jR d 1 

This formula does not cover two cases: d — 2,j — and d — 1,7 < 1/2. The estimate 
(118p in the borderline case d — 1,7 = 1/2 with a sharp constant was obtained in [TT] . 
Similar estimate for Jacobi matrices, which include the lattice one- dimensional Laplacian 
particular case, was proved in [12J. 
Our approach based on the annihilation of the underlying recurrent process allows us 
to obtain estimates on ^(V) for d — 1,2. Here we provide only two simple results. Other 
statements and a discussion are in section 9. The advantage of our estimate in the case 
d = 2, 7 G [0, 1] is related not to the fact that 7 can be zero, but to the independence of 
the constant from 7, contrary to (jTBl where q i7 — > 00 as 7 — > 0. 
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Consider one- dimensional operator (pQ) with a bounded potential: V(x) < A < oo. 
Then for 7 < 1/2 we have 

/OO POO 1 — 

■00 Jo s > 1 

Here T(7) is the gamma-function. This is a Bargman-type estimate for S^V). We will 
provide also a refined Bargman-type estimate. The next estimate is not very sharp in 
some cases (in particular it is much worse than in [11] when 7 = 1/2), but it is valid for 
each 7 > 0: 

/oo 
\x\V 1+ \x)dx, d=l. (19) 
-00 

The following estimate will be obtained for the two-dimensional Schrodinger operators: 
if < V(x) < 1, then there are some constants a±, a 2 such that for each 7 6 [0, 1], 

S y (V) < ai + a 2 - ln 2 (2 + \x\)dx. 

Jr 2 [n v^) 

Its analog is also valid in the lattice case without the assumption of the boundedness of 
the potential. 

The paper is organized as follows. Section 2 is devoted to the one-dimensional dis- 
crete and continuous Schrodinger operators. Two-dimensional continuous and lattice 
Schrodinger operators will be studied in sections 3 and 4, respectively. Proofs of some 
important lemmas will be given in the Appendix. Estimates from below are given in sec- 
tion 5 showing the sharpness of the results in sections 2-4. Operators on general discrete 
graphs are considered in section 6. In particular, it is shown there that one can't expect 
the translation-invariant estimates for N (V) to be valid when the underlying Markov pro- 
cess is recurrent. Sections 7 and 8 concern the fractional powers of the lattice Laplacian 
and Bessel operators. Lieb-Thirring estimates are studied in section 9. 

The authors are grateful to O. Safronov for productive discussions and to B. Simon 
for useful critical remarks. 



2 One-dimensional operators 



Theorem 2.1. The Bargmann (TT3j] and refined Bargmann ( fl^D estimates hold for the 
operator 

d 2 

H = — — -V(x) in L 2 (R). 
dx 2 

Remark. It was mentioned in the introduction that estimate f lT3|) is well-known. 
Proof. Consider a rank one perturbation Hi of the operator H = — which is 
obtained by imposing the Dirtichlet boundary condition at x = : 



#i = -^, D Bl = Qr(R)f\{y:y(0) = 0}. 
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Then 



where 



Then 



Thus 



N (V; H ) < N (V; H x ) + 1 < -L / V(x) [ Pl {t, x, x)dtdx + 1 

C{<7) J R Jo_ 

V(x) 

d d 2 



e 4t e 4t 
Pi(t, x, y) = — j= 7==r, x,y >0. 



Pi{t,x,x) = 1 Z_ 4t . (20) 
V47rf 



Similarly, the kernel R^\x,y) of the resolvent (—Hi — A) 1 satisfies 

i?\ J (x,y) = ^-j= , R\l(x,x) = -\x\. 

Since j °° Pi(t, x, x)dt = — R+q, the latter relation together with ([HD imply ( TTH]) . In- 
equality (fl~4"l) immediately follows from ((7j) and (120]) since the formula for pi above implies 
that 

1 POO J g"3r" 

pi(t,x,x)dt = \x\F( - ), F( 7 ) = / - — dr, 

V^XjaT ,/ 7 V47TT 

and F( 7 ) < 1 for all r > 0; Fh) < f°° -^=dr = when 7 > 1. □ 

\ 1 1 — — ' V ' / — J<y 4r\/47rr 4 V /7T7 ' — 

Consider now the same operator on the one- dimensional lattice: 

Hijj(x) = -Aip - V(x)i}) = 2ip(x) - i)(x + 1) - ip(x - 1) - V(x)il)(x) in L 2 (Z). 

The general solution of the equation Atp — Xip = 0, A > 0, on the lattice Z has the form 
ip = Ciaf + C^af , where a± t 2 are the roots of the equation a 2 — (2 + X)a + 1 = 0. If 
a = 2+A+ ^ A2+ — , A > 0, is the biggest root, then the solution of the equation 

(A — X)Rf\x, y) — 5(x — y) 

must have the form R^\x,y) = ca~\ x ~ y \ where the constant c can be easily found from 
the equation. This leads to 

If Hi is the lattice Laplacian with the Dirichlet boundary condition at x = and 

POO 

R^(x,y) = - P i(t,x,y)e- Xt dt 
Jo 
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is the kernel of its resolvent, then R^\x, y) = R^\x, y) — R,f\x, —y), and 

_ _1— 2|as| 

rW(x,x) = — , A>0. 

A v ' ; 2-(2 + A)a' 

We note that a ~ 1 + and 2 - (2 + X)a 2^X as A ^ +0. Hence, -R^ (x, x) = \x\ 

and therefore (Q, (fTTj) are proved for the one-dimensional lattice operator. 

In order to obtain a refined Bargmann estimate in the lattice case, we note that 

Pi(t,x,y)=p (t,x,y)-p (t,x,-y), where p (t, x, y) = ^ f e -«(i-«« *M*-v)*d<l>, 



Pi(t, x, x) = p (*, ^) - Po(*, x, -x) = p (t, 0, 0) - p {t, 2x, 0). 

The integral above can be expressed through the modified Bessel function. This allows 
to obtain the asymptotic behavior of po(t, x, 0) as t, \x\ —> oo. Another option is to apply 
Cramer's form of the central limiting theorem [9] (Ch. 16, 7) which leads to the following 
result: if t — )■ oo then 

p Q (t,x,0) = ^^(l + 0(-)), for \x\<t 2 '\ 

V Ant t 

\po\ < e~ ctl/ \ \x\ > t 2 '\ 

These formulas allow us to obtain the same estimate for pi(t, x, x)dx as in the contin- 
uous case, which leads to 

N Q {V)<d{a) J2 \x\V(x) + C 2 (<t) Yl x 2 vHx) + 1. 

^W>3 *:V(a)<^ 



3 Two-dimensional Schrodinger operator, continuous 
case 

This section is devoted to the estimate of Nq(V) for the two-dimensional operator ([1]) in 
L 2 (R 2 ). The rank one perturbation does not work in this case. We consider a soft killing 
by a potential instead. Let Hi = —A + q(x), x G R 2 , be the perturbation of the operator 
H = —A, where q(x) = 1 when \x\ < 1, q(x) = when \x\ > 1. Let p\ = pi(t,x,y) be 
the solution of the corresponding parabolic problem 

p u = Apx - g(x)pi, t > 0, pi(0,a;,y) = 5„(a;). (21) 

We need to show that the Markov process with the generator Hi is transient and we need 
a sharp estimate of pi as t — > oo. 
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at 



Theorem 3.1. The following estimate holds 

\pi(t, x, x)\ < ^ g — — when t > 7(2;) = max(l, a|x| 2 In \x\) (22) 
tin t 

for some a > 0. 

Proof. We will provide here only a sketch of the proof. The rigorous arguments will 
be given in the Appendix, Lemmas I10.ltil0.3l The first of these lemmas states that 

\pi(t,x,y)\ < —J-, M,M<2, t>2. (23) 
tin t 

In order to justify this estimate we solve (12T1) using the Laplace transform and arrive 

Pl = - J^\x,y)e xt d\, x^y, (24) 
where R ^ (x, y) is the kernel of the resolvent 

< = (A-g(a:)-A)- 1 , 

and the contour V consists of the bisectors of the third and second quadrants of the 
A— plane with the direction on V such that ImA increases when a point moves along V. 
We show (in the proof of Lemma [10. ip that the kernel (x, y) is bounded at A = 0, and 

(x, y) = a(x, y) + ^b(x, y) + O(j^), A — >• 0. (25) 

Obviously, 

e xt d\ = for t > 0. (26) 

Further, replacing T be a contour 7 around the negative semi-axis in the A— plane we 
obtain that 

1 r 1 r°° 1 1 
- e *'dA = / -l-e A *dA = / [- ; }e- at d(T 



In A J 1 In A J In cr + 7U In a — 7ri 

e - at da= — + 0{ — 5-), t>2. (27) 



\n 2 a + n 2 ' tln 2 t t\n 3 t 

The last two relations together with (1231) and (1241) imply (J23]). 

Lemma 110.21 extends (123]) for arbitrary x if i is large enough. It states that 

\pi(t,x,y)\ < gwM|±M ; | y | < 3/2, > 2 if t > max(l, a|x| 2 In |a;|) (28) 
tin t 
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for some a > 0. Note that In (2 + \x\) appears in the estimate for pi, and the estimate is 
valid only if t ^> \x\ 2 . The main step in the proof is a comparison of p\ with the function 

x) + b(j)(t, x), where 



iP(t, x ) = I v(\,x)e xt d\, v = v(X,x) = K (^\ x \\ (20) 

In A 



Here 



K(ji) = K Q ( fi ) = '^H^\if i ), /i>0, 



is the modified Bessel function (it is proportional to the Hankel function of the pure 
imaginary argument). 

Function if) satisfies the heat equation when \x\ > 2, vanishes when t = 0, \x\ > 2, and 

^M=2 = 7r27 + (7-r3T)> *->°°> ( 3 °) 

tin t tin t 

since expansion (125]) holds for v when |x| = 2. Thus there exist constants A and r such 
that Aip > pi when \x\ — 2, t > r. We choose to be the solution of the heat equation 
for |se| > 2 with zero initial data and the boundary condition at \x\ =2 being zero for 
t > t and one for t < r. If b is large enough, then pi < Aip(t, x) + b<p(t, x) at \x\ = 2, and 
therefore p\ < Atp(t,x) + b<f)(t,x) at \x\ > 2. It remains to estimate functions ip and \<p\ 
for \x\ ^> 2. 

The asymptotic behavior of ip for large values of t depends on the behavior of the 
integrand in ( |29|) as A — > 0. If we use only the main terms of the expansion of the 
modified Bessel function at zero, we obtain that 

^ ~ / , \ U e xt dX, t oo. 

Jr lnA 

After that, fl26|) and fl27j) imply fl28l) . A similar estimate is valid for (p. 

The last lemma in the Appendix extends (128]) (with an extra logarithmic factor in the 
right-hand side of the inequality) to arbitrary x and y and completes the proof of the 
theorem. □ 

This sketch of the proof avoids some difficulties and it does not explain the reason 
to introduce the condition t 3> \x\ 2 . The rigorous proof will be given in the Appendix. 
A similar result in the lattice case will be proved in the next section using probabilistic 
ideas. 

Theorem 3.2. Estimate ( T73j) holds for two-dimensional Schrddinger operators (TJP in 
L 2 (R 2 ). 

Proof. The Birman-Schwinger principle implies that 

N {V) = N (V + q; H x ) < N (2V; H,) + iV (2g; H 1 ). 
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Let us show that No(2q; Hi) < 1. Indeed, operator Hi perturbed by the potential —2q 
coincides with the operator —A — q(x), i.e., one needs to show that the latter operator 
has at most one negative eigenvalue. In fact it has exactly one eigenvalue, but we need 
only the estimate from above. The eigenfunctions i/j n {%) of this operator can be found by 
separation of variables, i.e., they have the form 

= {oin cos(n0) + /3„ sin(n8)) f n (r) , n>0, 

where r = \x\, 6 = arctan(y/x) and (— A — q(x))f n = ^ n f n - The corresponding eigenvalues 
are \ n = n 2 +^ n . Since the operator —A is strictly positive and q(x) < 1, we have j n > — 1. 
Thus X n > when n > and A > —1. In fact, the operator —A — q(x) can not have 
positive eigenvalues, i.e., this operator has only spherically symmetrical eigenfunctions 
/ = fo( r ) an d the corresponding eigenvalues A G (—1,0]. 

Let us show that the eigenvalue problem (—A — q(x))fo = Xfo has at most one eigen- 
value with a spherically symmetrical eigenfunction. We write the problem in polar coor- 
dinates: 

1 f°° 
- y"{r) - -y'{r) - q(r)y(r) = Xy(r), AG (-1,0], y(0) < oo, / ry 2 dr < oo. (31) 
r Jo 

The equation does not have nonzero solutions in L 2 (R + ) if A = 0. If y satisfies (l3Tj) with 
-1 < A < 0, then 



y = C l J (VlT\r), r<l; y = C 2 K (y/\X\r), r>l 



where Jo is the Bessel function and Kq is the modified Bessel function (Hankel function 
of the purely imaginary argument). Since Kq(g) ^ when a > and Jq{&) ^ when 
1 > cr > 0, the function y does not vanish at r > 0. If there exist two nonzero solutions 
of problem (|3T!) with different values of A, then at least one of them must change its sign 
on R + . Hence, N (2q; Hi) < 1, and the first inequality in ( !T5|) is proved. 

In order to prove the second inequality in (fT5|) we rewrite the first inequality in the 
form 

p /*oo p poo 

Nq(V) < c(cr) / V(x) / px(t, x, x)dtdx + c(a) / / V(x)pi(t, x, x)dtdx + 1, 

where Di = {x : > ^(x)}, D 2 = R 2 \Di and 7(2) is defined in ffT5l) or fl22l) with a — 1. 
Theorem 13.11 allows us to rewrite the first term in the right-hand side of the inequality 
above as the first term in (|15p . Hence, in order to prove the second inequality in (|T5|) it 
remains to show that 

(2 + \x\)V 

Pi(t, x, x)dt < C In ^ , x G D 2 . (32) 

a 
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We split the interval of integration in two parts (y^y, 1 + \J(1 + j(x), oo). We 
estimate function \pi | from above on the first interval by 1/ 4irt and apply Theorem I3.ll 
on the second interval. This implies that the left-hand side in fl32|) does not exceed 

. , a Cln 2 (2 + |x|) 
ln(l + 7) - In —— + ' ' 



V{x) m(l+7) 
The latter value can be estimated by the right-hand side in (l32|) since 

d ln(l + 7) < ln(2 + < C 2 ln(l + 7). 



□ 



4 A two-dimensional lattice operator 

We consider two-dimensional lattice operators 

Hip(x) = -Aip - V(x)ip in L 2 (Z 2 ) 

in this section. First we will prove (flUj) which justifies the Bargmann estimate (1121) 
(a = 0). Then we will prove the refined estimate (a > 0). 

The following two facts are the key starting points for the proof of f fTOj) . Let R^ (x, y), 
(x, y) be the kernels of the resolvents of the operators A = — Hq and —Hi, respectively, 
where —H\ is obtained from — Hq by imposing the Dirichlet boundary condition at the 
origin (annihilation of the Markov process at this point). Obviously, 

POO 

(x, y) = - pi (t, x, y)e- xt dt, A > 0, 



where p\ is the transition probability for the Markov process with the generator —H\. 
Since p\ > 0, one can pass to the limit as A — > in the relation above. Thus, the transience 
of the process x(t) for the operator Hi is equivalent to the condition \Rq (x, x)\ < 00, 
and it is sufficient to prove estimate (fTUj) for — R^\x, x). Secondly, the following relation 
is valid for the latter function: 

itf (*, x) = 2 A hm[i?f (0, 0) - Rf (x, 0)] . (33) 

Indeed, the kernel R\(x,y), A > 0, must have the form R\(x,y) = R\(x,y) + ci?°(x,0), 
where c = c(y) can be found from the condition R\(0,y) = 0. This immediately implies 



R^\x,x) = [Rf\x,x)- [I f^ ^ }, A>0, (34) 

R\ ( U ) U J 
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and 

I$\x,x) = lim [Rf\x,x) - [ ^fc°J ]. (35) 
A ^ +u it\ (0,0) 



Note that fj35 jl holds for general discrete operators. 

Formula ( 135]) for general discrete operators can be written in the form ( 133]) when the 
following three conditions hold: operator Hq is translation invariant, the Markov process 
with the generator — H is recurrent, and the Markov process with the generator —Hi is 
transient (in particular, if H Q = —A on Z 2 ). Indeed, ( 13"4"]) can be rewritten in the form 

(1) _ B(Rf\0,0) + Rf\x,0)) 
a(,) " i?f(0,0) 

where B = R ( °\o, 0) - R { °\x, 0). Since R ( °\o, 0) and Rf\x, 0) have the same sign (they 
are negative), the ratio satisfies 



Rf\0,0) 



(R ( °\0,0) + R { °\x,0)) 



€ [0,1]. 



Hence B is bounded when A — » 0. From here and lim^+o \R^\o, 0)| = oo, \Rq(x, x)\ < 
oo it follows that 

4% = lim [<(°.Q)];-[<^.o)] 2 = Um {B|2 _ = 21im 

i.e., 033]) holds. 

Theorem 4.1. Lei f/i = —A 6e i/ie negative Laplacian in L 2 (Z 2 ) with the Dirichlet 
boundary condition at x = 0. Then relations fTW . ZioZd. 

Proof. We will show that for each fixed x G Z 2 , 
R^\x,0) = — lnA + u(ar) + o(l) as A -> +0, where |«(a;)| < Cln(2 + |ac|). (36) 

This and (|33j) imply that R^\x,x) = 2[it(0) — u(x)\, and therefore \Rq\x,x)\ < Cln(2 + 
\x\). The latter is equivalent to (llOp and justifies (1121) . Thus it remains only to prove 

dsn]). 

The Fourier method applied to the equation (A — \)ip = 8(x) leads, for A > 0, to 



Rf\x,0) 



(2tt) 2 y [7ri7r]2 2 cos 0i + 2 cos 2 - 4 - A 

e i(*,<t>)d(f) 



(2tt) 2 7 h „ ]a 4 sin 2 & + 4 sin 2 f + A 



2 1 ikJ " i 2 
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where = (0i, 2 ) € [— 7r, 7t] 2 C -R 2 . We put here 

4 sin 2 ^ + 4 sin 2 ^ = |0| 2 + /*(0), |/i(0)| < C|0| 4 . 

The difference between (|37|) and the same integral with /i(0) = is 

1 f h((t))e l{x ^d(t) 



(2tt) 2 y hW)W]a [4 sin 2 f + 4 sin 2 § + A] [|0| 2 + A] ' 

The latter integral converges to a bounded function 

1 r h{(j))e i{x ^d(t) 
V[X) ~ J^f J { _^ ]2 4[sin 2 ^ + sin 2 f]|0| 2 

as A — >• +0, i.e., 

(2vr) 2 y hWi7r] 2 |0| 2 + A 
Obviously, the function 



Vi(x) = lim / 



[-7T,^\{H>1} 



l 2 + A 



is also bounded in x. Hence 



B§\ x ,0) = f ^JI+ W ( X )+ (1), A->+0, \w(x)\<C. (38) 
l Z7r J </|</>|<i \<P\ + A 

/" d0 



(2vr 

Note that 



_ ) + A -7r[ln(l + A)-lnA]. 



■/|*|<i 

Thus (138]) implies ( 136]) if it is shown that 

\F(x)\ < Cln(2 + |ar|), (39) 

where 

The function F, considered for all x E R 2 , depends only on r = \x\. One can replace x in 
the formula above by x = (r, 0). Then after passing to the polar coordinates a = |0|, 6 = 
arctan 02/0i, we obtain 



■|<i 191 Jo Jo 



This justifies ( 139]) and completes the proof of the theorem. □ 
The next statement provides the lattice analog of inequality (fT5|) . 
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Theorem 4.2. The following estimate is valid for the number of negative eigenvalues of 
the lattice two-dimensional Schrddinger operators: 

N (V)<d(a) i^ln 2 (2 + M) + C 2 (a) £ Vln^ + l, (40) 

x:V(x)<-f^r V ( x ) x-.VMy* 

where 7(2) = max(l, |x| 2 In \x\). 

Remarks. 1. The main difference between (|40p and (IT5|) is that the second integrand 
in (fT5|) contains V under the logarithm sign which is absent in ([40]) . Its presence in ( TT5|) 
is due to the non-integrability of the transition probability po(t,x,x) at t = for the 
Laplacian in i? 2 . 

2. After (1401) is proved, one can get a better estimate: 

M>flO<Ci(a) X) i^ln^+MHCMa) £ yin^ + iV+l, (41) 

where iV = : V(a;) > 1}. Indeed, let us introduce the potential V(x) which coincides 
with V at the points x where V(x) < 1, and V(x) = elsewhere. The operators H with 
the potentials V and V differ by an operator of rank N, and the difference between the 
numbers of their eigenvalues can be at most N. Thus estimate fT40l) for the potential V 
implies ( l4Tj) . 

Proof. We only need to show that 

Pl(t,X,x) < C^ 11 \ X \ wnen l > ry( X ), (42) 

tin t 

where pi(t,x,y) is the transition probability of the Markov process generated by A with 
the annihilation at the origin x = : pi(t, x,x) = P x {x(t) = x,x(s) ^ 0, s G [0,t]}. The 
estimate ( 14TJ]) follows immediately from fT42|) and (J7J). One can use the same approach to 
prove (T42l) as in the case of operators in R 2 , using the representation of p\ through the 
inverse Laplace transform of the resolvent (x, x) followed by the asymptotic analysis 
the corresponding integral. We decided to describe another approach to justify f T42|) using 
probabilistic ideas. However, we will provide only a sketch of the proof. 
We use formula fl34l) for the resolvent 



RV{x,x) 



POO 

/ e~ xt p 1 (t,x,x)dt. 
Jo 



The ratio gbjj^ m this formula has the following important interpretation. Let r 



min{t : x(t) = 0}. Then 

Rf\x,0) 



£f (0,0) 



POO 

E x e~ Xr = / e- Xs q x (s)ds, 
Jo 



(43) 
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where q x (s) is the distribution density of the random variable r (if the process starts at 
x.) The inverse Laplace transform of (I34p leads to the identity 

pi(t,x,x) = p (t,x,x) - / q x (s)p (t - s,x,Q)ds. (44) 



o 



One can derive from (1371) that 



R^(x,0) ~ — ln(V\\x\) as x G Z\{0}; i?i 0) (0, 0) ~ In V% 

Z7T 27T 

when \/A|:r| — >■ +0, A — > +0, respectively, i.e., 

>(o) 



— TTvi ~ , VAX ->■ +0. 

4°) (0,0) \ny/X 



(The latter formula has an analog in the continuous case if the annihilation occurs at the 
unit circle, and r is the time needed to reach the circle. Then the Laplace transform of 
the distribution density of r equals ^~t^^ ~ ^^ ). The latter formula and ( T43l) after 
the rescaling A = Ai|x| _a , a > 2, imply 

\ r a — 2 

E x e ^ — >■ as x — >■ oo. 

a 

Since the right-hand side above does not depend on Ai, it follows that n; converges in 
law to zero or infinity with the probabilities ^— ^ and -, respectively: 



, r , a; — 2 , t ,2 

|x| a a |x| Q a 

The latter relation after the substitution s = |x| Q , a > 2, (i.e., a = jjn^r) leads to 

2 In I x I s 

P x {r > s] ~ — when \x\ oo, -. — —> oo. (45) 

Ins \x\ 2 

After formal differentiation (in fact, we do not need to use the differential form of the 
relation above), we get (compare to Lemma fl 0.2 1) 



2 In |x| 
sin 2 s 



qx{Sj ~ ; 2 , S > 7(x), |x| -> OO. 



The main estimate (112]) follows from ( )45l) . formula ( 1441) and standard Gaussian estimate 
for po : 



„2 



Po(*,a;,0) = — (1 + 0-j=), \x\ > t^ 2 . 



□ 
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5 Estimates from below 



The goal of this section is to show that the estimates (llip - (fT5|l are sharp in the following 
sense: the operator has infinitely many negative eigenvalues in the case of any potential 
which decays at infinity a little slower (by a logarithmic factor) than in those estimates. 
To be more exact, the following theorem holds 

Theorem 5.1. Let H = —A — V(x) be a one- dimensional Schrddinger operator in L 2 (Z) 
or L 2 (R) with the potential V such that for some e > 0, 

x\ „ , , 

oo, 





x\ 




ln 1+£ ( 


1+1 


x\) 



— tj— ; — — V(x) = oo, or / — t—- ■ — — V(x)dx 

, ln 1+e (l + |x|) 1 J J R ln 1+E 1+ x) K } 



respectively. Then H has infinitely many negative eigenvalues (N (V) = oo). 

Let H = —A — V{x) be a two-dimensional Schrddinger operator in L 2 (Z 2 ) and 

£V(ar)=oo. ( 46 ) 

Then N (V) = oo. 

Proof. We will prove the first statement (d = 1) only in the lattice case. The 
continuous case can be treated similarly (and in fact, it is simpler). Consider sets I = 
h = {x : 2 k < x < 2 k+1 } C Z. Let 



a k 

Since J^a^ = oo, there exists an infinite sequence of values of k = kj, j = 1,2,..., for 
which 

a k > k- {l+£/2 \ k = k r (48) 

By taking a subsequence, if needed, we can guarantee that kj + i — kj > 2. Let L k = {x : 
2&-1 < x < 2 k+2 }, k = kj, be the union of l kj and two neighboring sets l k - The sets {L kj } 
do not have common points (except, perhaps, the end points). The first statement of the 
theorem will be proved if, for infinitely many sets L = L kj , we construct functions ip = ipj 
with the support in L and such that (Hif),ip) < 0. 
We will take 

ip = sin[— (x — a)}, x G L, ip = 0, x f. L, 

where \L\ = 2 k+2 — 2 k ~ 1 , k = kj, is the length of the interval between the end points of L 
and a = 2 kj ~ 1 is the left end point of the set L. The function ip is a sine function whose 
half-period is L and which is zero outside L. The l 2 (Z)— norm of this function for large 
L has order \/L/2: 

= ^4/2(l + o(l)), \L\^oo. 
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One can easily show that —A sin ax = a{a) sin ax, x G Z, where a is arbitrary and 
a(a) = 2 — 2 cos a ~ a 2 as a — > 0. Hence, 

7T 7T 
-A^ = 0"( prr)^ - Sin — (<5 (z) - <y b (x)), 

where is the delta function at the point y, and a, b are the left and right end points 
of L, respectively Thus, (—Aip,iJ)) = 1 1"0| | 2 , and therefore 

vr 2 

(-Ai/>,i/>) ~ -— \L\^oo. (49) 

Let us evaluate now 

(Vip,^) = J2v(x)ip 2 (x) > V(x)^ 2 (x). 

XdL Xdl 

Since ij+fe^ < C2 k k~ 1 ~ £ on l k and > 0, from gZD and (ggj it follows that 

> C x 2- k k £/2 > CilLl-Mn^ 2 |L|, |L| -> oo. 

Furthermore, is located far enough from the end points of and there exists c > 
such that if)(x) > c, i6 4. Hence, 



Cln £ / 2 


L\ 









Together with (149 p . this proves that (Hip,ip) < for large enough L. 

The proof of the one-dimensional statement of Theorem 14.11 is complete. 

Let us prove the statement of the theorem concerning the two-dimensional operators. 
As in the previous case, we will construct a sequence of functions ip = ipj(x), x G Z 2 , 
with non-intersecting finite supports and such that (Hip, ip) < 0. The functions ipj will be 
defined by induction as the restrictions of some functions — <Pj on the Euclidian space 
R 2 onto Z 2 C R 2 . Denote by Qk squares in R 2 for which \xi\, \x 2 \ < k. Let us define 
= 4>j +i while assuming that functions fy, j < j , have been constructed. We choose 
k large enough so that the supports of all the functions if)j already defined are located 
strictly inside Qk- We take k — 1 to define the first function 0i. The function = 4>j +i 
will be supported by a square layer P = Q 2 \ \ Qk with some I ^> k chosen below. Thus 
each layer Pj is split naturally in two parts, the interior part = Qi\ Qk and the 
exterior part P^ = Q 2 i \ Qi- We put = outside of P and = 1 on the interior part 
of the layer P. Then we split the exterior part P^ into four trapezoidal regions using 
diagonals of the square Q21 and define to be such a linear function in each of these 
trapezoidal regions that = 1 on the boundary dQi of the square Qi and = on dQ 2 i. 
Note that = on dP. 
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Let us estimate (— Aip,ip), x G Z 2 , from above. We will use notation dQk for the 
boundary of the square Qk C R 2 , and q for the union of the boundaries of the trapezoidal 
regions in P C R 2 constructed above. Since — Au = for any linear function u on Z 2 , 
the support of the function — Aip belongs to the set dQk [jdQk+i [jq, i.e., 

k-a^voi< E i A ^i' 

dQk U 9Q k+1 (J q 

since < ip < 1. Furthermore, \i^(xx) ~~ ^(^l)! — 1 fo r each pair of neighboring points 
xi,X2 G Z 2 , and therefore (see <Q) \ — Aip\ < 4, x G Qfc U dQ^+i- In fact, the latter 
estimate holds with 2 instead of 4, but we do not need this improvement. A better 
estimates holds on q. Since |V0| < 1/Z, we have | — Atp\ < 4/7, x G q. Taking into 
account that \dQk\ + \dQk\ < C\k and |Z| < C2I, we arrive at 

|(-A^,V)| < 4ciA; + 4c 2 . 

Note that the latter estimate does not depend on I. 

Obviously, (Vip,ip) > Ylx&pw V( x )- Assumption ( 1461) allows us to choose I such that 
the right-hand side of the latter inequality exceeds 4c\k + 4c 2 . Then (H^,ifj) < and the 
proof is complete. □ 

6 General discrete Schrodinger operators with recur- 
rent underlying Markov processes 

This section is devoted to a Bergmann type estimate for general lattice operators. We 
also will show here that shift-invariant estimates of the form (fT6l) can not be valid for 
operators with recurrent underlying Markov processes. 

Let X be a countable set and let Hq be a symmetric non-negative operator on L 2 (X) 
with matrix elements h(x,y), i.e. 

h ^(x) = E K x > y)i>{y), h ( x , y) = %, x )- 

It is assumed that 

h(x,y)<0 if x 7^ y, 2/) = 0; h(x,x)<co for all x G X. 

Obviously, operator H Q can be written in the form 

y£X:y^x 

The first two conditions above guarantee the existence and uniqueness of the Markov 
process x(t) with the generator — H and that the operator H is non-negative, and the 
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last one is needed to avoid a pathological behavior of the Markov process x(t). We also 
assume connectivity, i.e., X can not be split in two disjoint non-empty sets X 1 U X 2 in 
such a way that h(xi,X2) = for each Xi G Xi, x 2 G X 2 . If Hq = —A on Z d , we have 
h(x, x) = 2d, h(x, y) = —1 when \x — y\ — 1, h(x, y) — when \x — y\ > 1. 

Let po(t,x,y) be the transition probability, i.e., po is the kernel of the Markov semi- 
group e~ tH °, and let 

/•oo 

Rf\x,y) = - / p (t,x,y)e- xt dt (50) 
./o 

be the kernel of the resolvent = (—H — A) -1 of the operator —H . The connectivity 
assumption implies that po(t,x,y) > and R\ < for all values of the arguments. Since 
—H < 0, the function Rf\x,y) is analytic in A ^ (— oo,0]. We assume that the process 
x(t) is recurrent, i.e., 

roc 

/ Po(t, x, x)dt — oo. (51) 



The latter relation implies that 



\\m^Rf(x,x) = -oo. (52) 



The following result is a Bargmann type estimate for the lattice operator H = H 
V(x). Let us fix an arbitrary point Xq G X. Denote 



■ [flfWo)] 2 P (Q), 
^+o L 4 0) (x ,a;o) 



R = R( x , Xo ) = hmj ^ ' 0, { - R^(x,x)} 



and 

R = R(x,x ) = 2 lim [Rf\x,x ) - R { °\x ,x )} (53) 

A — ^~h0 

if the operator H is translation-invariant. 

Theorem 6.1. Let the Markov process x(t) with the generator —H be recurrent. Then 

1 ) the function R is finite for all x, Xo G X and positive for x ^ Xq (it vanishes if 

X = Xq), 

2) the following two estimates hold 



N (V) <Y,V(x)R(x,x ) + l, (54) 
N (V) < #{x G X : V(x) > 1} + V(x)R(x,x ) + l. (55) 

x:V(x)<l 

Proof. Let us denote by x(t) the Markov process x(t) with the additional condition 
of annihilation at the point x$ G X. From the connectivity assumption it follows that 
the conditional process x(t) with the annihilation at x G X is transient. With this fact 
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taken into account, the above formulas for R differ from ( )33|) . ( 133)) only by sign and by 
Xo playing the role of the origin x = 0. Thus, 

POO 

R(x, xq) = — R^\x, x) = / pi(t,x,x)dt, 

Jo 

where pi(t,x,y) is the transition probability for the process x(t). This implies the first 
statement of the theorem. The latter relation and ([7]) with a = lead to (134|) . In 
order to obtain (1331) . one can consider the potential V which is obtained by reducing 
the values of V to zero at all points where V > 1. Then (1331) follows from ( 13^1) for the 
potential V since the operators with these two potentials differ by an operator of the rank 
N = #{x G X : V(x) > 1}. □ 
Theorem 16.31 below shows that a space invariant estimate of N (V) can not be valid 
for the discrete operator H, but first we need to prove the following preliminary result. 

Theorem 6.2. For each y G X and a > 0, the operator H = H — a5 y (x) has a unique 
simple negative eigenvalue A = \(y,a) < 0. 

Proof. The uniqueness is due to the fact that if is a rank one perturbation of H Q . 
Let us show the existence of the eigenvalue. First we note that J2 x Po(t,x,y) = 1, and 
therefore (J50l) implies that 

^i?i 0) (x,y) = ^, A>0. (56) 

X 

Formula ( 130]) implies also that Rf\x,y) < 0. Thus from (1361) it follows that \R^\x, y)\ < 
j- for each x, y G X and A > 0. This and (1361) leads to the estimate ^2 x [Rf\x, y)] 2 < 
i.e., Rf{x,y)eL\X), y G X. 

We look for an eigenfunction in the form ip\(x) = R^\x,y), A > 0. Since (— H — 
A)^a = S y (x), will be an eigenfunction of H = H — a5 y (x) with the eigenvalue —A if 
—aR ® (y,y) = 1. The latter equation has a solution A = X(y, a) > for each y G X and 
a > due to (1321) and the relation lim^oo R^\y, y) = 0. □ 

The next theorem shows that for each 7 one can find a potential V > such that 
Ylxex ^(x) is arbitrary small and the operator H has infinitely many negative eigen- 
values. Hence, estimate (1T6|) can not be valid for the operator H. The potential V will 
be constructed when a uniformity condition on the unperturbed operator H holds. We 
assume that there exists an integer- valued metric d(x,y) on X (for example, I 1 — metric 
on Z d ) such that the following two relations hold: 

a) \R^\x,x)\ > /3(A) for A > and some /3(A) > 0, and /3(A) — > 00 as A — > +0; 

b) E x-.d(x,y)>r \ x i v)\ tends to zero uniformly in y when A > and r — > 00. 
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Theorem 6.3. Let conditions a),b) hold. Then for any sequence a n — > +0, one can find 
a set of points {x n G X} such that the operator 

oo 

H = H -J2aJ Xn (x) (57) 

n=l 

has infinitely many negative eigenvalues. 

Proof. In order to prove the theorem, it is sufficient to construct a sequence of 
compactly supported functions {ipk( x )} with disjoint supports such that 

(HMx),Mx)) <0. (58) 
For fixed y G X, a > 0, consider a "test" operator 

H = H(y, a) = H — aS y (x). 

Due to the previous theorem, this operator has a negative eigenvalue — Xo(y,a), where 
A = \o(y,a) > is the root of the equation —aR^\y,y) = 1. The corresponding 
eigenfunction can be chosen as 

VE,[<W)1 2 

Note that condition a) implies that Ao > Ao(a) > 0, where the lower bound Xo(a) does 
not depend on y. 

In order to complete the proof of the theorem, we will need the following lemma. 
Lemma 6.4. There exists a function r = r(a) such that the inequality 

(H$(x),$(x))<0, H = H(y,a), (59) 
holds for the truncated eigenfunction 

The important part of the statement of this lemma is that r is y— independent. The 
statement follows from the uniformity assumption. Indeed, 

(Hij>(x),ij>(x)) = -Xo(y,a) < -A (a) < 0. 

Hence, it is enough to show that 

\(Hif>,i/>) - (H$,$)\ = \(Hi/>,i/>-$) + (H(il> < when r>r(a). 
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Since the operator H is bounded in l 2 (X), \\ip\\ = 1, \\ip\\ < 1, it remains to prove that 
— — > uniformly in y when a is fixed and r — > oo (all the norms here and below 
are in l 2 (X)). It was shown in the proof of Theorem 16.21 that \R ®\x,y)\ < ~ for each 
x,y G X and A > 0. Thus, from condition b) and the estimate \o(y,a) > Ao(ct) > it 
follows that 

11(1 -I d (,,y)<r)R { S(x,y)f<l- J2 ftl/)h0 

x:d(x,y)>r 

uniformly in y when r — >■ oo. Condition a) implies that 

^[< ) (x,y)] 2 >[< ) (y,y)] 2 >/3(A )>0. 

a; 

This completes the proof of the lemma since 

II / 7l|2 _ IK 1 - / ^,2/)<r)^i° ) (a ; 



E,[<W)P 

Let us complete the proof of the theorem. We fix a\, calculate r = r(ai), select 
an arbitrary point X\ and chose ip\(x) to be the truncated eigenfunction of the "test" 
operator H(xi,ai). Other points x n , n > 1, will be chosen outside of the support of 
ipi. Thus inequality f )59|) with the "test" operator H(xx,cei) implies the same inequality 
for operator ( 157|) . i.e., (J5S]1 holds for tpi. Now we fix a 2 , calculate r = r(a 2 ), select 
an arbitrary point x 2 such that d(x 2 ,Xx) > r(a 2 ) + r (ai), and chose to be the 

truncated eigenfunction of the "test" operator H(x 2 ,ac2). The supports of functions ipi 
and ip2 are disjoint, and other points x n , n > 2, will be chosen outside of the supports 
of ipi,ip2- Hence, ( )58|) holds for %p 2 - The point x% will be chosen in such a way that 
d(x 3 ,Xi) > r(a 3 ) +r(«i) and d(x 3 ,x 2 ) > r(a 3 ) + r(a 2 ), etc.. This procedure allows us to 
construct the desired sequence {^(x)}. □ 



7 Fractional power of the lattice operator 

This section provides an illustration of the results on general discrete operators obtained 
in the previous section. It concerns an important class of non-local random walks x a (t) 
on the one- dimensional lattice. The one- dimensional lattice Laplacian 

-H = Ai/)(x) = ip(x + 1) + if>{x - 1) - 2ip(x) 

on l 2 (Z) is the generator of the symmetric random walk x(t) with continuous time. Let 
P t = e tA be the corresponding semigroup, P t ip = Y^yEzPi^^iV^iv)- ^ ne °P era tor A in 
the Fourier space L 2 [—tt,7i] acts as the multiplication by the symbol 

A(0) = 2(cos0- 1) = -4 sin 2 ^, G 
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Similarly, 

roo 1 

P t (0) = e - w i P A (0) = - / e - xt P t dt = - - , 2<j) , G [-7r,7r]. 

The main object that we study in this section is the fractional degrees Hq, a > 0, of 
the operator Hq = —A. After the Fourier transform, the operator Hq , its semigroup and 
the resolvent are the operators of multiplication by the symbols 



-A?=(4sin 2 £f, iV = e-< W r, R 



-1 



A + (4sin 2 |)« 

The following result is well-known in probability theory. 

Lemma 7.1. The operator —Hq = — (— A) a , a > 0, is the generator of a Markov process 
x a (t) on Z iff < a < 1. 

Proof. One needs only to prove the positivity of the kernel p a (t, x, y) of the semigroup 
P t)Q for < a < 1 and non-positivity of the kernel for a > 1. If < a < 1, then there 
exists [9] (Ch. 13, 6) the probability density g a ,i{s) > 0, < s < oo, (which is called the 
stable law with the parameters a and (3 = 1) such that 

POO 

e- x " = / e' xt g ajl (t)dt. 
Jo 

The second parameter /? in the two-parametric family of the densities g a ^ characterizes 
the symmetry of the density. If (3 = then g a ,o( s ) = g a ,o(~ s ), if /5 = 1, < a < 1, then 
9a,i(s) = 0, s < 0, g a>1 {s) > 0, s > 0. 

Using the probability density g a ^i, one can represent Pt, a in the form 

poo poo 

Pt, a = e~ t{ ~ A)a = / e tl/asA g aA (s)ds = / P^ s g a>1 {s)ds, 
Jo Jo 

i.e., the kernels p a and p of the operators P tta , P t are related by 

POO 

p a (t,x,y)= / p(t 1/a s } x } y)g a ,i(s)ds. 



This implies the positivity of p a . 

In order to show that the semigroup P t)a is not positive when a > 1, we note that 
the function h(t,(f)) = Pf iQ ,(0) has the following property: h"(0) = 0. Hence, its inverse 
Fourier transform h(t, z) = f_ h(t, <j))e lzr ^d<f) satisfies ^2 zeZ z 2 h(t, z) = 0, which shows 
that p a = h(t, x — y) can not be non-negative. □ 

Lemma 7.2. For each a G (0, 1] and t — > oo, 

r(l/(2a)) 

p a [t,x,x) ~ tl/(2a) , c Q 2 ^ . 
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Corollary. The random walk x a (t) is transient for < a < 1/2 and recurrent for 
1/2 < a < 1. The formula above indicates that p a (t,x,x) has the same asymptotic 
behavior as the transition probability for "nearest neighbors random walks" (defined by 
the standard Laplacian) when the dimension d equals I /a. A similarity between the 
long range 1-D ferromagnetic interaction and high-dimensional local interaction (similar 
to noted above) was discovered by Dyson j8]. This similarity was a foundation for the 
introduction of the hierarchical lattice. We will discuss the spectral properties of the 
hierarchical Dyson's Laplacian elsewhere. 

Proof. This statement is a simple consequence of the Laplace method applied to the 
integral 

pJt, x,x) = — r e-< 4sin2 ~ — f e-M^dh t -> oo. 

2tt J-n 2tt J_ n 

□ 

Remark. Similar calculations give a more general result. If x a (t) is a random walk 
on Z with the generator — (— A)°, then 

> 4> in law as t — v oo, (60) 

where (p has a stable distribution g2a,o(s), s G R, (symmetric stable law with parameters 
2a, (3 = and characteristic function (Fourier transform) equal to e~ x Indeed, 

E e * 1/a = e V ?7^y ->• e , t — >■ oo. 

Formula f )60|) means that, after rescalling, the lattice operator (— A) a approximates the 
fractional power of the continuous Laplacian, i.e., the random walk x' a (s) = x a (st) jt x l a 
(which is the rescalling of x a (s)) converges to the stable process x* a (s) on R with the 
generator -(-^) a . 

The following theorem is an immediate consequence of the standard CLR estimate 
(JJ]) with a = 0, Lemma 17.21 and finite rank perturbation arguments (see Remark 2 after 
Theorem I4.2[) 

Theorem 7.3. (Transient case) Consider the Hamiltonian 

H a = -(-A) a - V(x) on l 2 (Z), V(x) > 0, 

with < a < 1/2 (i.e., the dimension d — - > 2). Then there is a constant C = C(a) 
such that 

N (V) < #{x e Z : V(x) > 1} + C(a) ^ V^(x). 

x:V(x)<l 

Remark. The constant C(a) can be evaluated. One can show that C(a) = 0( x \ ) 
as a — > 1/2. 
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Let us consider now the recurrent case: a > 1/2. First we calculate the regularized 
resolvent (|53|) : 

i r 1 - e i<t>x 



Ro(x,0) = lim [Rx, a (x, 0) - ii! A)a (0, 0)] = lim / 

A->+0 A^+0 Z7T 



a^+o 2tt A + (4sin 2 f) Q 

= |im 2 r shf(fe) 2 r 

A^+OTrio A + (4sin 2 §)« Wo (4sin 2 |)« 

A simple analysis provides the following asymptotics of the regularized resolvent when 
|x| — )■ oo : if a > 1/2, then 

~ 4 f°° sin 2 z 1 

R (x,0) ~ c a \x\ 2a -\ c a = - -^—dz = 0{- -) as a ^ 1/2. 

If a = 1/2 (a borderline case which corresponds to a Cauchy-type random walk with the 
generator —(—A) 1 / 2 ), then 

Rq(x, 0) ~ — In 1x1. 

7T 

Hence, Theorem 16.11 implies 

Theorem 7.4. (Recurrent process) There exist constants C a ,C such that 

N (V)<#{xeR: V(x)>l} + C a ^ Vix)^ 2 ^ 1 + 1, ^ < a < 1, 



x:V(x)<l 



N (V) < #{xG R : V(x) > 1} + C V(x)\n(2 + \x\) + l, a = 1/2. 

x:V{x)<l 

8 Bessel operators 

This section concerns another class of one-dimensional operators which may have an 
arbitrary positive spectral dimension. These operators Bd on the half line R + are defined 
by the radial part of the Laplacian: 

d 2 d-ld 

d = TH + T> r > °- 

dr A r dr 

We consider arbitrary (not necessarily integer) d > 0. The operator can be represented 
in the form 

B^^I—^T^ d> ^ ^ 
r a 1 r ar 

i.e., it is symmetric in L 2 ({0,oc),r d ~ 1 dr). 
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If d > 2, then the operator 5^ in £ 2 ([0, oo), r d ~ 1 dr) is self-adjoint (we do not need to 
impose a boundary condition at r = 0). The diffusion process bd{t) with the generator 
Bd is also well-defined since the point r = is not accessible from any initial point r > 0. 
If d < 2, the situation is different. The equation B^ip = has two bounded linearly 
independent solutions, ipi — 1)^2 

= r 2 ~ d , i.e., due to the Weil criterion, we have the 
limit circle case near r = 0, and a boundary condition is needed to define a self-adjoint 
operator. One can impose the Dirichlet boundary condition at r = which corresponds 
to annihilation of the process bd(t) at r = 0. Another option is an analog of the classical 
Neumann boundary condition: lim^+o r ^'(r) = (see [13]). 

We will consider the Schrodinger operator Hd = —Bd — V(r) with the Dirichlet bound- 
ary condition at r = if d < 2 and without a boundary condition if d > 2. The process 
bd(t) in both cases is transient, and our main concern is to obtain an exact formula for 
the transition probability pd(t,a,r) with respect to measure dfi = r d ~ l dr. Denote by I v 
the modified Bessel function of order v. 

Lemma 8.1. The process bd{t) has the following transition density if d > 2 or d < 2 and 
the Neumann boundary condition at r = is imposed: 

p d (t,a,r) = (2t)- 1 e-^(ar) 1 - d / 2 / (i/2 _ 1 (^). 

If d < 2 and the Dirichlet boundary condition is imposed, then 

-d/2 



p d (t,a,r)=p°(t,a,r) = (2t)- 1 e"^ (ar^h-a^). 



The first formula can be found in [13J. The second formula can be proved similarly. 
From Lemma [8.11 it follows that 

A—d 

Pd(t,r,r) ~ t oo, d > 2, p° (t, r, r) ~ Cd^—^ t ->• oo, d < 2. 

Applying the CLR estimate, we obtain 
Theorem 8.2. Ifd>2, then 

POO 

N {V) < c{d) / V d/2 r d - l dr. 
Jo 

If d < 2 and the Dirichlet boundary condition at r = is imposed, then 

N?(y)< Cl (<T) f V(r)r 2 - d dr + c 2 (a) f V 2 - d/2 {r)r 4 - 2d dr. 

Jr:r 2 V>a Jr:r 2 V<a 

The standard rank one perturbation arguments imply that the last estimate with 
constant one added to right-hand side is valid for N^(V). 
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9 Lieb-Thirring sums 



The results of this section are based on two known formulas for the Lieb-Thirring sums for 
the general Schrodinger operators H = Hq — V(x) on a complete a-compact metric space 
X. The first formula is valid under the same assumptions, that are needed for formula 
OH) (in particular, the transience of the underlying Markov process is required), and has 
the form 

S 7 {V) < [ V 1+1 {x) [ p (t,x,x)dtfx(dx). (62) 



cur 



x 



V(x) 

d 



Note that formula (fTBl) for the operator H = —A — V(x) in R , d > 3, is an immediate 
consequence of Q62]) . The second formula is valid under the same conditions, but the 
transience requirment is replaced by the convergence of the following integral: 

t~~ f po(t,x,x)dt < oo. 



If the latter integral converges, then 



■V-.m < 2 ^j7 / V{x) I r^p (t,x,x)dtfi(dx), (63) 



where IX7) is the Gamma-function. Note that fl6"3"|) implies ffTSl) for the operator H = 
—A — V(x) in R d when | + 7 > 1, i.e., the case d = 1, 7 < 1/2 is still not covered by 
(!62|) .Q63l). Estimates for this case will be obtained below. Our approach allows us also 
to obtain new estimates in the cases when (|62j) or (|63|) hold, and in some cases these 
new estimates are better than ( |62l) or (1631) . For example, our estimates on S^(V) for the 
two-dimensional Schrodinger operators are uniform in 7 e [0, 1], while the right-hand side 
in fl63|) goes to infinity for these operators when 7—7-0. 

While estimate fl62|) can be found in many papers starting from the original paper 
by Lieb and Thirring [18] (see also [22], [25], [IS]), we didn't find a reference for (1631) 
(similar formulas can be found in [22]). Thus we decided to give a brief proof of it. Let 
Ne(V) = #{Ai <-E,E> 0}. We have 



SytY) = 7 



poo POO 

/ E 1 ~ 1 N E (V)dE < 2 7 / E^- l Ti(V[(H + E)' 1 - (H + V + E)~ l ])dE 
Jo Jo 

/"OO /"OO 

< 2 7 / E^TriV / [ e -^ Ho+E ^ - e-^ Ho+v+E ^}dt)dE 
Jo Jo 

poo 

= 2 7 r( 7 ) / t-^Tr(V[e- tHo - e-^ Ho+v ^])dt. 







The right-hand sides here and in fl63|) coincide, and this justifies (1631) . 

Consider now a Schrodinger operator H = Hq — V(x) on a metric space X such that 
the Markov process x(t), generated by — H , is recurrent, and a point xq is accessible 
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from any initial point. Then the process Xi(t) with annihilation at the moment of the 
first arrival to Xq is transient. For example, H Q can be a negative lattice Laplacian on 
Z d , d < 2, the general discrete operator discussed in section 6, or the generator of a 1-D 
diffusion process, say, H = — x G R, or H = —B^d < 2, see (l6Tl) . Let H\ be the 
generator of the process x\(t). It is given by Hq with the Dirichlet boundary condition at 
x : ip(x ) = 0. Let pi(t,x,y) be the transition probability for the process Xi(t). 

We will assume additionally that the potential is bounded: V(x) < A. This implies 
that the ground state A (U) is bounded from below, A (U) > —A. Since the operator 
H = H\ — V(x) is a rank one perturbation of H = Hq — V(x), the eigenvalues of the 
operators H and H alternate. Hence, the bound for the ground state and estimates 
([62]) . (|63|) applied to the operator H lead to the following statement. 

Theorem 9.1. Let Sy(V) be the Lieb-Thirring sum for the Schrddinger operator H = 
Hq — V(x), where H is an operator which satisfies the conditions described above, and 
V(x) < A. Then the following two estimates hold 

SJV) < A 7 + -^- / V 1+7 (x) / Pl (t,x,x)dtfi(dx), (64) 
c{a) J x J« 

V(x) 

S 7 (V) < A 7 + 2 ^_y / V(x) I r"'p l {t,x,x)dtfi{dx), (65) 



Remark. The second formula can be applied to the potentials which decay at infinity 
slower than in (164"]) . but it worsens when 7 — > 70 and t~ l0 pi is not integrable at zero. 

Let us apply the latter theorem to the one-dimensional Schrddinger operator on R. 
We choose xo = 0. Using formula ( IS"4"1) we arrive at (compare with ( JT4l) ) 

Theorem 9.2. Let H = — — V(x) be the one- dimensional Schrddinger operator on 
L 2 (R) and < V(x) < A. Then 

S y (V) < A 7 + -^-[/ \x\V l+ \x)dx + [ x 2 V 3/2+ ^(x)dx}. (66) 

C \ a ) Jx 2 V(x)>a V " 71 " Jx 2 V(x)<u 

Using f )65|) and the same arguments as in section 2 we obtain the following statement. 



Theorem 9.3. Let H = —rns — V{x) be the one- dimensional Schrddinger operat 
L 2 (R) and < V(x) < A. Then for any 7 < 1/2, 

S^V) < A 7 + -^r[ci / \x\ 1 ~ 2l V{x)dx + c 2 [ x 2 V 3/2+ ^(x)dx], 

C \ a ) Jx 2 V(x)>a Jx 2 V(x)<a 



or on 



where 



-1 



c - 7r(7) r i ~ e ~- ,/ , 



n J s^)/ 2 ' A (1 + 27)^+^' 
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Note that c\ — > oo when 7 — > 1/2. 

Let us turn now to the two-dimensional Schrodinger operator on L 2 (R 2 ). We can 
not use Theorem 19.41 in this case since each point xq G R 2 is not accessible for the 
two-dimensional Brownian motion. Formula ( 1631) provides an estimate for S^(V) with a 
constant which blows up when 7—7-0. One can obtain a better estimate for small 7 using 
annihilation due to a compactly supported potential g(x) introduced in section 3: q — 1 
for \x\ < 1, q = for \x\ > 1. The main theorem is a consequence of fl63l) . where the 
operator iJ = —A — V(x) is considered as the perturbation of —A + q by the potential 

-y-g. 

We will assume that the potential V is bounded (otherwise the formula is too cumber- 
some). Then one can use the scaling and reduce the problem to the case when V(x) < 1. 

Theorem 9.4. Let H = —A — V(x) in L 2 (R 2 ) and < V(x) < 1. Then there exist 
constants ax, a 2 < 00 such that for each 7 G [0, 1], 



Sj(V) < ai + a 2 / V +1 [ X ^ ln 2 (2 + \x\)dx. 

JR 2 ln T777T 



lR 2 m V ( x ) 

Proof. Let Hi = — A + q(x) and let pi(t,x,y) satisfies 
dp 1 

— + ifipi = 0, t > 0, pi(0,x,y) = Sy(x). 

Then 

ln 2 (2 + |x|) . . 1 

Pi(t, x, x) < c — 7- ' u , t > 2, x,x)<— < t < 2. 

cm r 47rr 

Indeed, the second estimate is due to the fact that p\ < po, and it is valid for all the 

values of the arguments. The first estimate was proved in Theorem 13.11 (see also Lemma 

I10.3j) for t > |x| 2 ln(2 + \x\). It remains to note that the validity of the first estimate for 

2 < t < |x| 2 ln(2 + |x|) follows immediately from the fact that ln pjj^ < C < 00 for those 

t. 

We apply formula fl63l) with a = 4 for the operator —A + q(x) perturbed by the 
potential — V — q and arrive at 

S y (V) < ^7T^ I (V(x) + q(x)) [ ri pi {t,x,x)dtdx. 

R? 



c(4) 

Since 



V(x)+g(x) 



V(x) + q(x) 



we have 

/ x , Z" 00 ln 2 (2 + Ixl) , 

t 7 Pl(t, X, XjCtt < Cn / n dt 

v ; _ J A t^n^t 
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c ln 2 (2 + \x\) f°° dt c ln 2 (2 + \x\) 



A~i J A tln 2 t A" 1 In A 

Thus 

SJV) < C x / ln 2 (2 + \x\)G(V(x) + q(x))dx, where G(z) = -r—r, 0<z<2. (67) 
! R 2 In * 



The following subadditive inequality holds for the function G: there exists a constant C 2 
such that 

G{z x + z 2 ) < C 2 (G( Zl ) + G(z 2 )) for all Zl , z 2 > 0, Zl + z 2 < 2, < 7 < 1. 

The last two inequalities imply the statement of the theorem. □ 
Theorem 19.41 remains valid in the lattice case without the assumption of the bound- 
edness of the potential. 



10 Appendix 

Here we obtain estimates on the solution pi of problem ff2TT) as t ^> \x\ 2 — > 00 which 
provide a rigorous proof of Theorem 13.11 

Lemma 10.1. The following estimate holds j or p^: 

\pi{t,x,y)\ < -Sr-, |x|,|y|<2, t > 2. (68) 
tin t 

Proof. Consider the operator A — q(x) in L 2 (R 2 ). This operator is negative, and its 
spectrum coincides with the semi-axis (— 00, 0], i.e., the resolvent = (A — q(x) — A) -1 
is analytic in A G C = C \ (—00, 0]. We apply the Laplace transform to (I2T1) and arrive 
at 

pi = - / i?i 1} (x, y)e A 'dA, a > 0, x ^ y, 

J a—ioo 

where R^\x,y) is the kernel of operator R^\ i.e., 

(A-q(x)-\)R^\x,y) = 6 y (x). (69) 
We deform the contour of integration in the integral above: 

Pl = - ^R^(x,y)e xt d\, x^y, (70) 

where the contour T consists of the bisectors of the third and second quadrants of the 
A— plane with the direction on T such that ImA increases when a point moves along 
T. This deformation of the contour is possible since for each fixed x ^ y the function 
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R^\x,y) is analytic in A G C , decays as |ImA| — > oo and is bounded as A G C, A — > 0. 
The boundedness will be justified below when the asymptotic behavior of R^\x,y) as 
A ->■ is established (see fITT]) . (JT2])). 

Our next goal is to find the asymptotic expansion of R^ (x,y) as A — > 0. Let us 
represent R^ (x, y) in the form 

R<k\x, y) = ^\n\x-y\+ u(X, x, y), \y\ < 2, (71) 

where x 6 Co°; X = 1 f° r M < 3, X = for |a;| > 4. Since R^\x,y) has a logarithmic 
singularity when a: — y y, the function u is bounded when \y\ < 3 and A > is fixed. 
We intend to show that 



1 

hTA'" v ' : hr A' 



u = u (x,y) + J:rT u 1 (x,y) + rTT v(\,x,y)), \u \, |tti|, \v \ < C, (72) 



when \x\, \y\ < 2, A G C , |A| < 1. 
We put (ITTj) in (1691) and arrive at 

(A - g(a;) - X)u = }\ + Xf 2 , 

where 

1 1 ^(x) 

fi = -Vx • V ln \x - y\ + tH a x ~ 90*0) ln l x ~ 2/1' h = — z~ ln \x-y\- 

7T 27T 27T 

Obviously ||/i|| L 2, H/2IU2 < C, \y\ < 2, and f 1 = f 2 = for \x\ > 4. 

In order to describe the behavior of the function u when |A| < 1, we consider the 
truncated resolvent 

= TzR^T, : L 2 com (R 2 ) -> L 2 oc (R 2 ), XeC' = C\ (-00, 0], 

of the operator A-q(x). Here T x : L 2 com (R 2 ) L 2 (R 2 ) and T 2 : L 2 (i? 2 ) L 2 oc (i? 2 ) are the 
imbedding operators. Thus, the truncated resolvent is defined on the space of the square 
integrable functions with compact supports, and the images are restricted to bounded 
regions in R 2 . The following facts can be found in [52"j [55] : the truncated resolvent R^ 



A 

is analytic in C, admits a meromorphic continuation on the Riemannian surface of the 
function ln A and has the following asymptotic behavior at the origin 

^ = Ao+ ]^X Al + lr7A 5(A) ' A(h Au B{X) : L ^ R2 ) L L(R 2 ), (73) 

where the operators A , A 1: B(X) are bounded and ||5(A)|| < C as A — > 0, A G C. The 
constant C in the latter estimate depends on the size of the supports of the functions / 
in the domain of the operator B and the size of the domain in R 2 where the functions Bf 
are considered. 
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The validity of (173"]) needs an explanation. In fact, ( 173]) holds in our case only be- 
cause q(x) > 0. In the case of more general potentials q(x) (or more general operators), 
expansion ( 173]) has a more complicated structure [33J involving half-integer powers of A 
and polynomials of In A. Formula (1731) is an immediate consequence of this more general 
result in the case when \\R^\\ < C as A 0, A G C". The uniqueness of the solution of 
the problem 

(A-q(x))w = 0, x G R 2 ; \w\ < C, \Vw\ < C\x\~ 2 as x -» oo (74) 

implies (see [33]) the boundedness of \\R^ as A — > 0, A G C", and therefore it leads to 
( 173]) . It remains to note that 

=< (A - q(x))w, w >— - / (\Vw\ 2 + q(x)\w\ 2 )dx 

Jb? 

for solutions of (174")) . Thus w = if g(x) > and q is not identically equal to zero. Hence, 
(J73D holds. 

We imply the standard a priory estimates (for the Sobolev space Hf oc (R 2 )) followed by 
the Sobolev imbedding theorem and replace the space Lf oc in ( 173"]) by the space C(\x\ < 2) 
of continuous functions on the disk \x\ < 2. This and the formula u = R^\fi + Xf 2 ) 
complete the proof of expansion (1721) . 

We will also need an estimate on the function u in (17T1) when ImA — > oo. Since 
u = R x (fi +A/2) and the norm of the resolvent R^ does not exceed the inverse distance 
from the spectrum, we obtain that 

\\u\\ L *m <C, \y\<2, A G Ti = T f|{A : |A| > 1}. 

This and the standard a priory estimate for elliptic equations imply that ||w||ir 2 (ii a ) < C\X\, 
where H 2 is the Sobolev space. Using the Sobolev imbedding theorem, we arrive at 

|u| < C|A|, \x\, \y\ < 2, A G r x . (75) 

Let us substitute expression £2D for R^\x,y) in ([70]). Since J r e xt dX = for t > 0, 
we arrive at 

Pi = — J u(X,x,y)e xt dX, x^y. 

Using ( J75l) . we obtain 

Pi = - [ u(X,x,y)e xt dX + 0(e- £t ), \x\, \y\ < 2, x + y, t>2, 

where r 2 = T \ Fx = T f]{X : |A| < 1} and the estimate of the remainder is uniform in |x| 
and \y\ (one can take e = y/2/2). The integral J r2 e xt dX can be evaluated. Since it has 
order 0(e~ et ), t — > 00, expansion ( 1721) implies 

\pi\<C x \( T^e xt dX\+C 2 [ -L^ e ReXt d\X\+C 3 e- £t , \x\,\y\<2, x^y, t>2. 
Jr 2 ln A Jr 2 I In A| 
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One can replace here the integrals over T2 by the same integrals over T plus the terms 
that decay exponentially as t — > 00. Then we make the substitution A = [i/t in the second 
integral which leads to the estimate of this integral by 7^7, t > 2. Thus 

M <Ci| [ -^-e Xt dX\ + —^-, \x\,\y\<2, x^y, t>2. 
Jr hi A t In t 

This formula together with (1271) prove (1681) for x ^ y. The latter restriction can be 
dropped since pi is continuous when t > 0. □ 
The next lemma provides an estimate on pi when |x| is not bounded. 

Lemma 10.2. Let \y\ < 3/2, \x\ > 2 and t > a\x\ 2 In |x| /or some a > 0. T/ien 

\ Pl (t,x,y)\<^^, C = C{a). (76) 
rln r 

Proof. Consider the following function v = v (A, X ) = ^M where 

^)=ir (^) = y^ 1) M, /i>0, 

is the modified Bessel function ( it is proportional to the Hankel function of the purely 
imaginary argument). The function v is the exponentially decaying, as |x| — > 00, solution 
of the problem 

(A - \)v = 0, Id > 2, v\ lxl=2 = h(X) := K( f^\ \eC = C\ (-00, 0]. 

In A 

Let 

if)(t,x) = / v(X,x)e M d\, 



where T is the contour introduced in the proof of the previous lemma. Since v is analytic 
in A £ C and decays exponentially when | argA| < 37r/4, |A| — > 00, the integral J r vd\ 
vanishes, and ip is the solution of the problem 

V>t = AV>, |ar|>2; ^\ ]x]=2 = / h(X)e xt dX; V|t=o = 0. 



It will be shown below that the following estimates are valid for the function ip: 

^lM= 2 = ^ + 0(^), t^oo, (77) 

\ib\ < C ^ n \ x \ Ixl > 2, t > alxl 2 In Ixl. (78) 

In particular, f J77|) and Lemma 110.11 imply the existence of constants A and r such that 
A0 > pi when |x| = 2, t > r. 
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We will also need the solution of the following parabolic problem 

4 = A0, |z|>2; 0| N=2 = <?(*); 0| t =o = 0, (79) 
where g(t) = for t > r, g(t) = 1 for t < r. We chose a constant b large enough, so that 

b > max (pi + 

|z|=2, t<T 

Then from the maximum principle it follows that Aip + b<j) > p\ for all (t, x) : t > 0, |se| > 2. 
It will be also shown that estimate (1751) holds for the function when t > t: 



|0| < — ■ — 2~~"> 1^1 > 2, t > max(r, a|x| 2 In (80) 
tin t 

Since is bounded in any bounded region, ( 1501) implies the same estimate without the 
restriction t > r. The latter estimate together with ( 175j) imply ( 176|) . Thus the proof of 
Lemma [10.21 will be complete as soon as ( 177|) . (175]) . and ( 150]) are justified. 

Let us justify (177]) . From the logarithmic behavior of the function K(\/X) at zero and 
the exponential decay at infinity it follows that 

K(V\) = -\nV\ + a + f{\), |/(A)| < C|Aln A|, A G T, (81) 

where a is a real constant. We took into account here (and below) that In A ^ when 
A G T. Thus 

h(X) = = - 1 + hTX + g(A) ' |g(A)l < C ' |A| ' A e r " 

Expansion (1771) is a consequence of (127]) and the following two obvious relations: 

e xt d\ = for t > 0, | / g(A)e A 'dA| < — . 

Jr t 

Let us justify ( ]75j) . Expansion ( 15T]) implies 

= /■ g(yXN) At = r i "-Hxi At /■ « e , di (82) 

r 7 r In A 7 r 2 In A ; 7 r In A v ; 

Inequality (175]) with arbitrary t > 2 holds for the first integral in the right-hand side 
above due to (I2"7j) . In order to estimate the second integral on the right, we note that 

, /(AM 2 ) , ^ < C | A | N 2 + c , In |x| , 



In A In A In A 

Thus, the substitution A — > [i/t implies that the second integral does not exceed 



\x\ 



\x\ 2 In \x\ 
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Hence, flZE} holds if |x| > 2, ^ < cminf-^T, rM with some c > 0. It remains to note 
that the latter restrictions on x, t with some c = c(a) follow from those imposed in the 
statement of Lemma 110.21 since the equation t = a|x| 2 In |x| implies that lnt/ln|x| — > 2 
as | a; | —> oo. Hence, (j?8]) holds. 

In order to prove the validity of (15U1) we solve problem ( 17§j) using the Laplace transform. 
This leads to 

' <fi{X,x)e xt dX, 4>{X,x) - 



The contour of integration above can be replaced by T when t > t : 



2tt 



4>(\,x)e xt d\, t > r. 



r 



bmce AT(2V^) ^ in C", expansion flHTl) leads to the following representation of <f)(\,x) 
when A G T 2 = T f){\ : |A| < 1} : 

J _ - MVAW) + a + /(ALP) + = t + 2rlnW + C(A) ,„ w + Q{ / W) , 



-lnv / A + a + 0(AlnA) In A -2a In A 

The estimate 



^, f(A|a;| 2 ) . . Clnlxl . . 

0( - \ ' r e Af rfA < x > 2, t > a x 2 ln b 

In A ; 1 ~ tln 2 t — i i i 



f 2 



was proved in the process of evaluating of the second term in the right hand side of (1821) . 
Further, 

|ln|x| [ 0{X)e xt dX\<C\n\x\ [ \X\e ReXt d\X\ < C \n\x\ [ |A|e ReA *rf|A| = C lu '' r| 
and the integral J T e xt dX can be evaluated, which implies that 

f e xt dX = 0(e- £t ), t^oo. 



The last three estimates prove the validity of (18011 for the function 

1 



Zix 



(j)(X,x)e xt dX 
r 2 



if we take into account that 



1 e xt dX = [ — ^ e xt dX - [ e xt dX, r x = T\r 2 , 



r 2 lnA — 2a J r \iaX — 2a J ri \nX — 2a 
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where the following estimates hold for the terms in the right-hand side: 

I / — e xt d\\ <C I e KeXt d\\\ = 0(e- £t ), t ->■ oo, 

In A -2a 1 " J Tl 

— e xt d\\ < — t > 2. 

In A -2a 1 ~ tlnV 

The latter inequality can be proved absolutely similarly to (1271) . Hence, it remains to 
show that flHUl) holds for the function 



k 1 



2ti 



(fr(\,x)e xt d\, t>r. 



We note that K(2y/\) does not vanish on T and decays exponentially at infinity. This 



implies that | ^wtP I < C when A G r 1; \x\ > 2. Since ReA < —y/2/2 on r 1; it follows 



that 

|^| < C I \ l ~^ XT e xt d\\ < C{T)e- tV ~ 2 ' 2 t > r. 
Jr 1 A 

□ 

Lemma 10.3. Let \x\ > 2 and t > a|x| 2 ln \x\ for some a > 0. Then 

\ Pl (t, x ,x)\<^^-, C = C{a). (83) 
tux t 

Proof. While x can be an arbitrary point of R 2 below, it is always assumed that \y\ > 2 
(otherwise the estimate of p\ is provided by Lemma fl 0.1ft . Let \ — x{^i x ) e t > 0, 
be a cut-off function such that x = 1 when \x\, t < 1, x = when |x| > 3/2 or £ > 2. We 
are looking for pi(t,x,y) with |y| > 2 in the form 



e 4 * 

Pi(t, x, y) = p (t, x, y) - x{t, x)p (t, x, 0) + z(t, x, y), p (t, x, y) -- 



Ant 

Then z is the solution of the problem 

z t = Az- q(x)z - f, x G R 2 , z(0, x, y) = 8 y (x), (84) 

where 

/ = q(x) \x(t, x)po(t, x, 0) - po(t, x, y)\ + 2V X ■ Vp (*, a?, 0) + (Ax - Xt)Po(t, x, 0) G C°°. 
Note that the following estimate holds for the function h = x(i, x)po(*, 0)— po(t, x, y): 

C\x\ 2 

\h(t,x,x)\ < — ^- when |x| > 2, t > 2. 

v 
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In particular (see the proof of (175j) ). from here it follows that 

( a i In \cc\ 

\h(t,x,x)\ < 2 when \x\ > 2, t > a\x\ 2 In |x|. (85) 
tin t 

Hence, it is enough to prove the statement of the lemma for the function z instead of p\. 
We need some estimates on the function / in order to estimate z. Note that 

l„|2 i„ „,|2 



!= ^-^-h^' l> - 2 ' (86) 

Assuming that \y\ > 2 and taking into account that |x| < 1 on the support of q and that 
/ is bounded when \y\ > 2, t < 2, we obtain 

/ e « e « f 1 — e « A 1 — e « 

i/c*.-.y)i* < c+jf i-sr—sri*^ —s*—* 

The next estimate for / is valid when \y\ > 2, t > ||y| 2 ln|y|. From (156]) it follows 
that 2 

|/(t,x,y)|<^|L, |j/|>2, t>2, 
and therefore (see the arguments in the proof of (1781 ) 

l^ Mil ft | ,n , . . 

\f(t,x,y)\< \ 3 y >2, t>-y 2 lnj/. (88) 

tin i z 

We solve problem (1541) using the Duhamel principle and arrive at 

\z(t,x,y)\< / / Pi{s,x,x )\f(t- s,x ,y)\dx ds. 

JO J\x \<3/2 



Hence ^ 

\z(t,x,x)\< / / pi(s,a;,a;o)|/(t - s,x ,x)\dx ds 

JO J\x \<3/2 

+ / / Pi(s,x,x )\f(t — s,Xo,x)\dx ds, \x\ > 2. (89) 

Function pi(s,x,xo) can be estimated through p for all values of s, and therefore 

C 

\pi(s,x,x )\ < when \x — xq\ > 1/2. 

1 + s 

This and (155!) imply that the first term in the right-hand side of (159]) does not exceed 
c ^ when |x| > 2, t > ^|-y| 2 ln \y\. The same estimate is valid for the second term 

t In t 1 ' ~ ' — 2 1 " i 1 " 1 

due to (1571) and (ITS"]) (where one also needs to keep in mind that pi(s,x,y) = px(s,y,x)). 
Hence ([55]) is proved for z. □ 
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